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Hydrodynamics  of  Falling  Mine  in  Water  Column 

Peter  C.  Chu,  Anthony  F.  Gilles,  Chenwu  Fan,  Peter  Fleischer 


Abstract—  The  hydrodynamic  features  of  a  falling  mine  into 
the  water  column  is  investigated  experimentally.  The 
experiment  consisted  of  dropping  three  cylindrical  model 
mines  of  various  lengths  into  a  pool  where  the  trajectories 
were  filmed  from  two  angles.  The  controlled  parameters  were, 
mine  parameters  (length  to  diameter  ratio,  center  of  mass 
location),  and  initial  conditions  (initial  velocity,  and  drop 
angle).  Results  indicate  that  center  of  mass  position  has  the 
largest  influence  on  the  mines’  trajectory  and  that  accurate 
trajectory  modeling  requires  the  inclusion  of  both  momentum 
and  moment  equations.  A  statistical-dynamic  model  has  been 
established  to  predict  the  trajectories  of  the  falling  mines. 

Key  Words —  Mine  impact  burial,  momentum  of  mine, 
moment  of  momentum  of  mine,  mine  impact  burial  prediction 
model,  mine  drop  experiment  (MIDEX). 

1.  INTRODUCTION 

N  December  31,  1991  the  Union  of  Soviet  Socialist 
Republics  (USSR)  effectively  ceased  to  exist  under 
international  law  and  the  Cold  War  ended  (Fischer 
1999).  In  response,  the  Navy-Marine  Corps  team  developed 
a  new  strategic  concept,  "From  the  Sea”  (FTS),  which 
provides  a  framework  for  Naval  operations  into  the  21st 
century.  FTS  effectively  shifted  operational  focus  from  blue 
water  operations  to  sea-based  power  projection  into 
regional  littoral  areas.  FTS,  the  1994  revision  ’’Forward 
From  the  Sea”  (FFTS)  and  ’’Operational  Maneuver  from  the 
Sea”  (OMFTS)  all  provide  guiding  principles  for  sea-based 
power  projection  to  regional  littoral  areas  of  the  world. 

One  of  the  greatest  threats  to  U.S.  sea-based  power 
projection  in  littoral  areas  is  the  naval  mine.  Mines  were 
first  developed  in  1776  and  have  been  used  in  most  major 
conflicts  since.  Today,  an  estimated  50  countries  possess 
some  sort  of  mining  capability.  (Lehr  2000)  Mines  can  be 
used  in  both  offensive  and  defensive  roles.  Offensively, 
they  can  be  placed  in  enemy  waters  or  nearby  sea-lanes  in 
order  to  harass  military  and  commercial  shipping. 
Defensively,  they  can  be  used  to  delay  or  prevent 
amphibious  assaults  or  to  deny  command  of  the  sea.  The 
Wonsan  Korea  Mine  Crisis  and  Iraq’s  use  of  mines  during 
Desert  Storm  provide  excellent  examples  of  the  value  of  the 
naval  mine  as  a  defensive  weapon.  Shortly  after  the 
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October  1950  Wonsan,  Korea  mine  crisis,  then  Chief  of 
Naval  Operations  Admiral  Forest  Sherman  exclaimed, 
"when  you  can’t  go  where  you  want  to,  when  you  want  to, 
you  haven’t  got  command  of  the  sea.  Command  of  the  sea  is 
the  bedrock  for  all  of  our  war  plans.  We  have  always  been 
submarine-conscious  and  air-conscious.  We  have  now 
commenced  to  become  mine-conscious  beginning  last 
week.”  (Boorda  1999). 

Within  the  past  15  years  three  U.S.  ships,  the  USS 
Samuel  B.  Roberts  (FFG-58),  Tripoli  (LPH-10)  and 
Princeton  (CG-59)  have  fallen  victim  to  mines.  Total  ship 
damages  were  $125  million  while  the  mines  cost 
approximately  $30  thousand.  (Boorda  1999)  Mines  have 
evolved  over  the  years  from  the  dumb  "homed”  contact 
mines  that  damaged  the  Tripoli  and  Roberts  to  ones  that  are 
relatively  sophisticated  -  non-magnetic  materials,  irregular 
shapes,  anechoic  coatings,  multiple  sensors  and  ship  count 
routines.  Despite  their  increased  sophistication,  mines 
remain  inexpensive  and  are  relatively  easy  to  manufacture, 
upkeep  and  place.  As  such,  they  are  an  efficient,  yet  potent, 
force  multiplier  and  are  widely  available  to  any  country  or 
group  who  has  a  modest  ability  to  purchase  them. 

Naval  mines  are  characterized  by  three  factors:  position 
in  water  (bottom,  moored,  rising,  floating),  method  of 
delivery  (aircraft,  surface,  subsurface)  and  method  of 
actuation  (acoustic  and/or  magnetic  influence,  pressure, 
contact,  controlled).  The  littoral  battlespace  is  divided  into 
five  regions  based  upon  water  depth.  Within  each  of  these 
regions  naval  forces  can  encounter  multiple  types  of  threats 
(Fig.  1).  The  littoral  regions  are: 

•  Deep  Water  (DW).  Water  depths:  >300  ft.  Threat: 
mainly  moored  and  rising  mines,  although  a  few 
large  bottom  mines  exist. 

•  Shallow  Water  (SW).  Water  depths:  from  40  to 
300  ft.  Threat:  bottom,  moored  and  rising. 

•  Very  Shallow  Water  (VSW).  Water  depths:  from 
10  to  40  ft.  Threat:  bottom,  moored,  rising  and 
controlled. 

•  Surf  Zone  (SZ).  Water  depths:  <  10  ft.  to  the 
beach  itself.  Threat:  same  as  VSW  but  land  mines 
and  obstacles  can  also  be  encountered. 

•  Craft  Landing  Zone  (CLZ).  Water  depths:  the 
beach  itself.  Threat:  conventional  land  mines  and 
obstacles.  (U.S.  Naval  Mine  Warfare  Plan  2000) 

The  shift  in  focus  from  the  blue  water  to  the  littoral  has 
brought  many  new  challenges  to  the  warfighter.  The 
greatest  is  what  impact  will  the  highly  variable  littoral 
environment  have  on  future  operations,  particularly  mine 
countermeasures  (MCM).  The  most  influential 
environmental  parameter  to  successful  MCM  operations  is 
the  local  bathymetry  character  of  the  bottom.  This  key 
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parameter  often  determines  whether  an  area  should  be 
swept  or  hunted.  Bottom  clutter  in  the  form  of  rock 
outcrops,  coral  reefs,  man-made  debris  and  irregularities  in 
slope  provide  false  contacts  or  create  shadow  zones  that 
increase  overall  clearance  times.  Soft  bottom  sediments 
such  as  marine  clays  and  silts  cause  a  high  degree  of  mine 
burial  upon  impact.  These  buried  or  partially  buried  bottom 
mines  are  of  greatest  concern  to  the  MCM  planner. 

Environmental  data  collection  in  a  potential  adversaries' 
littoral  region  is  often  hampered  by  inaccessibility.  This 
lack  of  accurate  data  causes  a  certain  degree  of  uncertainty 
in  MCM  planning.  As  a  result,  several  numerical  models 
for  predicting  mine  impact  burial  (IB)  have  been 
developed.  The  most  promising  IB  model  was  originally 
developed  by  Amone  and  Bowen  in  1980.  Later 
improvements  by  Satkowiak  (1987),  Hurst  (1992),  Chu  et 
al.  (1998,  2000a,  b)  and  others  have  resulted  in  the  current 
version,  IMPACT  25.  IMPACT  25  creates  a  two- 
dimensional  time  history  of  a  cylindrical  mine  as  it  falls 
through  air,  water  and  sediment  phases.  IB  prediction  is 
largely  calculated  from  the  marine  sediment  characteristics 
and  mine  impact  orientation  and  velocity.  Chu  et  al.  (2001) 
found  overprediction  of  burial  depth  by  IMPACT25  after 
analyzing  the  synchronized  environmental  and  mine  burial 
data  collected  during  the  Mine  Impact  Burial  Prediction 
Experiment  (MIBEX)  in  the  Monterey  Bay  in  June  2000. 

The  essential  elements  of  the  mine  impact  burial  model 
translate  into  the  science  and  engineering  of  hydrodynamic 
process  of  a  falling  object  and  of  sediment  transport.  Any 
solid  object  falling  through  fluid  (air  and  water)  should 
obey  two  physical  principles:  (a)  momentum  balance,  and 
(b)  moment  balance.  The  current  IMPACT25  model  only 
considers  the  momentum  balance  of  the  mine  and 
disregards  the  moment  balance  of  the  mine.  Such  an 
incomplete  hydrodynamics  in  the  model  leads  to  unrealistic 
prediction  of  the  mine  falling  in  the  water  (no  helicoidal 
motion).  If  considering  momentum  and  moment  balance, 
the  falling  object  should  have  a  helicoidal  motion.  Without 
the  helicoidal  motion,  the  IMPACT25  may  over-predict  the 
impact  burial  depth. 


Figurel.  Littoral  mine  threat  (from  Rhodes,  1998). 


In  this  study,  a  nonlinear  dynamical  system  is  established 
for  the  movement  of  a  nonuniform  (center  of  gravity  not 
the  same  as  the  center  of  volume)  cylindrical  model  mine  in 
the  water  column.  A  mine-drop  experiment  was  conducted. 
The  experimental  results  show  the  nonlinear  characteristics 
of  the  trajectory  pattern.  The  data  collected  from  the 
experiment  can  be  used  for  mine  impact  burial  prediction 
model  development  and  verification. 

2.  DYNAMICS 

A.  Earth  Coordinate  System 

Two  coordinate  systems  are  used  to  describe  a  mine 
falling  through  the  water  column:  earth  and  relative 
coordinates.  The  earth  coordinate  system  is  fixed  to  the 
Earth  surface  with  horizontal  sides  as  x  and  y-axes  (along 
the  two  sides  of  the  pool),  and  vertical  direction  as  z-axis 
(upward  positive,  Fig.  2).  Suppose  a  mine  falling  into  the 
water  column.  The  mine  rotates  around  its  main  axis  (r{) 
with  an  angle  y/\  and  an  angular  velocity  of  Q .  Its  position 
is  represented  by  the  center  of  mass  (COM),  and  its 
orientation  is  represented  by  two  angles:  iff2  and  y/3  (Fig. 

3).  Here,  y/2  is  the  angle  between  the  rraxis  and  the 
horizontal  plane;  and  y/3  is  the  angle  between  the 
projection  of  the  main  axis  in  the  (x,  y)  plane  and  the  x-axis. 
The  angle,  y/2  +  n  /  2,  is  usually  called  attitude. 

The  relative  coordinate  is  rigidly  connected  with  the 
mine.  The  origin  (O)  of  the  relative  coordinate  system 
coincides  with  the  center  of  mass  (COM);  the  rraxis  is 
along  the  central  line  of  the  mine;  the  r2-axis  is 
perpendicular  to  the  plane  constructed  by  rraxis  and  z-axis 
(rrz  plane);  and  the  r3-axis  lies  in  the  (rrz)  plane  and  is 
perpendicular  to  rraxis.  The  selection  of  axes  (x,  y,  z)  and 
(ru  r2,  r3)  follows  the  right-hand  rule.  Let  V*  =  (Fi*,  F2*, 
F3*)  be  the  three  components  of  the  velocity  of  COM,  i.e., 
be  the  origin  velocity  of  the  coordinate  system  (r1?  r2,  r3). 
The  geometric  center  (GC)  is  located  at  (x>  0,  0).  For  GC 

below  COM,  x  >  0,  and  for  GC  above  COM,  x  <  0.  The 
relative  coordinate  system  (ru  r2,  r3)  is  obtained  through 
the  translation  and  two  rotations  y/2  and  y/3  of  the  earth 

coordinate  system.  Let  the  position  vector  P  be  represented 
by  PE  and  PB  in  the  earth  and  relative  coordinate  systems, 
PE  and  PB  are  connected  by 


cos  ys3 

-  sin  \f / 

o" 

cos  y/2 

0 

sin^2 

X 

P  = 

E 

sin^3 

COS  Iff 

0 

0 

1 

0 

P  + 

B 

ym 

0 

0 

1_ 

-  sin  \j/2 

0 

cos  iff  ^ 

z 

(i) 


2 


Journal  of  Counter-Ordnance  Technology  (Fifth  International  Symposium  on  Technology  and  Mine  Problem) 


3 


where  (xm*,  ym*,  zm*)  represent  the  position  of  COM  in  the 
earth  coordinate  system. 

B.  Nonlinear  Dynamical  Equations 

Any  solid  object  falling  through  a  fluid  (air  and  water) 
should  obey  two  physical  principles:  (1)  momentum 
balance  and  (2)  moment  of  momentum  balance.  Let  Vw  = 

(Vwu  Vwi,  vw3)  be  the  water  velocity,  and  ( co\,(o\,(0 *)  be 

the  components  of  the  angular  velocity,  referring  to  the 
direction  of  the  relative  coordinate  system.  The  independent 
and  dependent  variables  are  made  non-dimensional  by: 


where  g  is  the  gravitational  acceleration,  and  L  the  length  of 
the  mine.  The  non-dimensional  momentum  equations  for 
COM  are  given  by  (Mises  1959) 

dV  p-p  F 

— L  +  (oV3  -  0)V2  = - —  sin  i/si  +  — 1 — ,  (2a) 

dt  2  3  3  2  p  2  pgU 

dV  F* 

— 2-  +  ^  =— 2— ,  (2b) 

dt  3  1  pgU 

dV  p-p  F 

— -- coVx  -  - -cos^;  + — - — ,  (2c) 

dt  p  ^  pgU 

where  n  is  the  volume  of  the  mine,  p  the  mine  density, 

pw  the  water  density,  (Fx*,  F2 *,  F3*)  the  components  of 

water  drag.  The  non-dimensional  equations  of  the  moment 
of  momentum  for  axial  symmetric  mine  are 
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where  is  the  distance  between  COM  and  geometric 

center  (GC),  (Mi*,  M2*,  M3*)  the  components  of  the 
moment  due  to  water  drag,  (Ju  J2,  J3)  the  three  moments 
of  gyration, 

J,  =  { (r'  +  r3)dm*,  J2=  J (r32  +  r2)dm*, 
J3=\(r2+r2)dm*.  (4) 

The  orientation  of  the  mine  (y/2,y/3)  is  determined  by 
du /  du / 

— L  =  a)2 ,  cosi/f2 — -  =  0)3.  (5) 

dt  dt 
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Figure  2.  Earth  coordinate  system. 
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Figure  3.  Cylinder  orientation  and  relative  coordinate  system. 

The  eight  non-dimensional  nonlinear  equations  (2a-c),  (3a- 
c),  and  (5)  are  the  basic  system  for  determining  the  mine 
movement  in  the  water  column. 

3.  MODEL  MINES 

A.  Description 

Three  model  mines  were  used  for  the  drop  experiment  at 
the  Naval  Postgraduate  School  swimming  pool.  All  had  a 
circular  diameter  of  4  cm,  however  the  lengths  were  15,  12 
and  9  cm  respectively.  The  bodies  were  constructed  of  rigid 
plastic  with  aluminum-capped  ends.  Inside  each  was  a 
threaded  bolt,  running  lengthwise  across  the  mine,  and  an 
internal  weight  (Fig.  4).  The  internal  cylindrical  weight 
made  by  copper  was  used  to  vary  the  mine’s  center  of  mass 
and  could  be  adjusted  fore  or  aft.  The  center  of  gravity  of 
the  model  mine  is  the  origin  of  the  body  fixed  coordinate 
system. 

The  model  mine  is  composed  of  six  uniform  cylindrical 
parts  (Fig.  4):  (a)  a  plastic  hollow  cylinder  C(1)  with  mass  of 
mi,  outer  and  inner  radii  of  R\  and  R2 ,  length  of  (L  -2/i), 
and  the  center  of  gravity  for  the  part  (COMP)  to  be  at  its 
geometric  center  located  along  the  rraxis  at  r  =  x  \  (b)  an 

aluminum-capped  left  end  (Fig.  5)  solid  cylinder  C(2)  with 
mass  of  m2,  radius  of  Ru  length  of  lu  and  COMP  located 
along  the  rraxis  at  rx  =  LI 2-  l\!2  +j;  (c)  an  aluminum- 
capped  right  end  solid  cylinder  C(3)  with  mass  of  m3,  radius 
ofRu  length  of  /1?  and  COMP  located  along  the  rraxis  at  rx 
=  L/2- 1\!2  -  x  ;  (d)  a  cylindrical  thread  C(4)  with  mass  of  m4, 
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radius  of  R3,  length  of  ( L  -  2/i),  and  COMP  located  along 
the  rraxis  at  rx  =  x  \  (e)  a  cylindrical  threaded  bolt  C(5) 
with  mass  of  m5 ,  out  and  inner  radii  of  R2  and  R3,  length  of 
/2,  and  COMP  located  along  the  rraxis  at  rx  =  LH  -  x  -  h- 
l2l 2;  (f)  an  adjustable  copper  cylindrical  weight  C(6)  with 
mass  of  m6,  outer  and  inner  radii  of  R2  and  R3,  length  of  /3, 
and  COMP  located  along  the  rraxis  at  rl  =  S  +  x,  where 
S  is  the  distance  between  the  COMP  of  the  adjustable 
weight  and  the  geometric  center  of  the  model  mine. 


B.  Moments  of  Gyration 

Since  the  six  parts  (all  mines)  all  have  uniform  mass 
distribution,  the  moments  of  gyration  for  these  parts  are: 

m  wi,  9  9  m  m .  9  m ,  9 

+  K ),  j?  =  —  K ,  ■/’  =  — ^  > 
2^2  2 

rzn  AW,  9  ^  mA  9  .  9  9x 

J  }  =  —R%  J  =—R%  J[  =  —  (R2  +  ), 


y=y  =^[3r2i+3r22+(l-2iJ], 
12 


j<2»  =  j<2>  =  ^(3i?2  +  /2)  jl3>  =  j(3)  =  — 43£2  +  /2), 

2  3  v  1  1  2  3  v  1  1 

12  12 


Figure  4.  Internal  components  of  the  model  mine. 
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Figure  5.  Internal  structure  of  the  model  mine. 
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/s>  =  j<!»  =^(3i?J  +/2), 
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where  the  superscripts  for  the  moments  indicate  the 
cylindrical  parts.  The  resultant  moments  of  gyration  is 
computed  by 


J2=J,  =  +m{z2  - L~zY  +>w( - L  +  x)2 

/=!  "  2  2 


2  L  l  2 

+m4X  +m5{—  -  X~lx  )  +aw6(J  +  j). 
2  2 


(7) 


According  to  the  definition  of  COM,  the  coordinate  of  GC 
(  x )  is  determined  by 


\mfL/2-li-l2/2)-mS\ 

2>, 


(8) 


which  indicates  how  the  adjustable  weight  determines  the 
location  of  COM  for  the  model  mine. 


C.  Model  Parameters 
C.l.  Length/Diameter  and  Density  Ratios 

Our  goal  was  to  choose  a  scale  that  was  somewhat 
representative  of  the  real  world  ratio  of  water  depth  to  mine 
length,  but  at  the  same  time  would  be  large  enough  to  film 
and  would  not  damage  the  pool's  bottom.  The  model  mines 
were  based  on  the  realistic  assumption  that  a  3  m  mine  is 
laid  in  water  depths  of  45  m,  thus  producing  a  15:1  ratio. 
The  depth  of  the  pool  is  2.4  m.  From  this  ratio,  the  length 
(Z)  of  the  model  mine  is  chosen  as  1 5  cm.  The  addition  of  a 
12  and  9  cm  length  allowed  for  later  comparison  of  the 
sensitivity  of  water  phase  trajectory  to  the  ratio  of  mine 
length  over  diameter.  The  outer  radius  of  the  model  mine  is 
2  cm.  The  length/diameter  ratios  {LID)  are  15/4,  12/4,  and 
9/4.  The  corresponding  density  ratios  {pi p  )  are  1.70, 
1.68,  and  1.88,  respectively. 


C2.  x  -Value 

In  each  of  the  three  model  mines,  the  location  of  the 
weight  (i,e.,  the  value  of  S)  is  adjustable.  Use  of  (7) 
location  of  the  COM  {x  -value)  can  be  determined  (Table 
1).  The  positive  x  -value  indicates  that  COM  is  below  the 
geometric  center,  and  the  negative  x  -value  indicates  that 
COM  is  above  the  geometric  center. 
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Table  1.  Model  L/D,  density  ratio,  and  X -value  (unit:  cm). 


Mine 

L/D 

p!  P, 

Z\ 

*2 

Z0 

Z-i 

Z-2 

1 

15/4 

1.70 

1.85 

3.69 

0 

-1.85 

-3.69 

2 

12/4 

1.68 

1.21 

2.43 

0 

-1.21 

-2.43 

3 

9/4 

1.88 

0.68 

1.37 

0 

-0.68 

-1.37 

4.  MINE  DROP  EXPERIMENT 

An  mine  drop  experiment  (MIDEX)  was  conducted  at 
the  NPS  swim  pool  in  June  2001.  The  purpose  of  the 
experiment  is  to  collect  data  about  mine’s  motion  in  the 
water  column  for  various  combinations  of  the  model  mine 
parameters.  It  basically  consisted  of  dropping  each  of  three 
model  mines  into  the  water  where  each  drop  was  recorded 
underwater  from  two  viewpoints.  Figure  6  depicts  the 
overall  setup.  The  controlled  parameters  for  each  drop 
were:  LI D  ratio,  X  -value,  initial  velocity  (Vin),  and  drop 
angle.  The  Earth’s  coordinate  system  is  chosen  with  the 
origin  at  the  comer  of  the  swimming  pool  with  the  two 
sides  as  x-  and  y-axes  and  the  vertical  z-axis.  The  initial 
injection  of  mine  was  in  the  (y,  z)  plane  (Fig.  2). 

A.  Initial  Velocity 

Initial  velocity  (Vin)  was  calculated  by  using  the 
voltage  return  of  an  infrared  photo  detector  located  at  the 
base  of  the  mine  injector.  The  infrared  sensor  produced  a 
square  wave  pulse  when  no  light  was  detected  due  to 
blockage  caused  by  the  mine’s  passage.  The  length  of  the 
square  wave  pulse  was  converted  into  time  by  using  a 
universal  counter.  Dividing  the  mine’s  length  by  the 
universal  counter’s  time  yielded  Vin.  The  mines  were 
dropped  from  several  positions  within  the  injector 
mechanism  in  order  to  produce  a  range  of  Vin.  The  method 
used  to  determine  Vin  required  that  the  infrared  light  sensor 
be  located  above  the  water’s  surface.  This  distance  was  held 
fixed  throughout  the  experiment  at  1 0  cm. 

B.  Drop  Angle 

The  drop  angle  (initial  y/^ )  was  controlled  using  the 

drop  angle  device.  Five  screw  positions  marked  the  15°, 
30°,  45°,  60°,  and  75°.  The  drop  angles  were  determined 
from  the  lay  of  the  pool  walkway,  which  was  assumed  to  be 
parallel  to  the  water’s  surface.  A  range  of  drop  angles  was 
chosen  to  represent  the  various  entry  angles  that  air  and 
surface  laid  mines  exhibit.  This  range  produced  velocities 
whose  horizontal  and  vertical  components  varied  in 
magnitude.  This  allowed  for  comparison  of  mine  trajectory 
sensitivity  with  the  varying  velocity  components. 


Figure  6.  MIDEX  equipment. 


C.  Methodology 

For  each  run  the  mines  were  set  to  a  X  -value.  For 
positive  X  -value,  the  mines  were  placed  into  the  injector 
so  that  the  COM  was  located  below  the  geometric  center. 
For  negative  x  -value,  the  COM  was  located  above  the 
geometric  center  to  release.  A  series  of  drops  were  then 
conducted  in  order  of  decreasing  mine  length  for  each 
angle.  Table  2  indicates  number  of  drops  conducted  for 
different  drop  angles  and  x  -value  for  L/D  =  15/4.  Number 
of  drops  for  other  L/D  ratios  (12/4,  9/4)  is  comparable  to 
that  for  L/D  ratio  of  15/4.  All  together  there  were  712 
drops.  Each  video  camera  had  a  film  time  of  approximately 
one  hour.  At  the  end  of  the  day,  the  tapes  were  replayed  in 
order  to  determine  clarity  and  optimum  camera  position. 

Table  2.  Number  of  drops  conducted  for  different  drop 
angles  and  X  -values  for  L/D  =  15/4. 


wT 

15° 

o 

o 

45° 

o 

O 

75° 

z2 

13 

15 

15 

15 

12 

Z\ 

9 

15 

15 

15 

9 

Z0 

12 

14 

15 

18 

6 

Z- 1 

0 

6 

6 

6 

0 

Z-2 

2 

6 

6 

0 

0 

5.  DATA  RETRIEVAL  AND  ANALYSIS 

A.  Data  Retrieval 

Upon  completion  of  the  drop  phase,  the  video  from 
each  camera  was  converted  to  digital  format.  The  digital 
video  for  each  view  was  then  analyzed  frame  by  frame  (30 
Hz)  in  order  to  determine  the  mine’s  position  in  the  (x,  z) 
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and  (y ,  z)  planes.  The  mine's  top  and  bottom  positions  were 
input  into  a  MATLAB  generated  grid,  similar  to  the  ones 
within  the  pool.  The  first  point  to  impact  the  water  was 
always  plotted  first.  This  facilitated  tracking  of  the  initial 
entry  point  throughout  the  water  column.  The  cameras  were 
not  time  synchronized;  thus,  the  first  recorded  position 
corresponded  to  when  the  full  length  of  the  mine  was  in 
view. 

B.  Source  of  Errors 

There  were  several  sources  of  error  that  hindered  the 
determination  of  the  mine's  exact  position  within  the  water 
column.  Locations  above  or  below  the  camera's  focal  point 
were  subjected  to  parallax  distortion.  Placing  the  cameras 
as  far  back  as  possible,  while  still  being  able  to  resolve  the 
individual  grid  squares,  minimized  this  error.  Second,  the 
background  grids  were  located  behind  the  mine's  trajectory 
plane.  This  resulted  in  the  mine  appearing  larger  than 
normal.  This  error  was  minimized  by  not  allowing  the 
plotted  points  to  exceed  the  particular  mine's  length.  Third, 
an  object  injected  into  the  water  will  generate  an  air  cavity. 
This  air  cavity  can  greatly  affect  the  initial  motion, 
particularly  at  very  high  speeds  (hydro  ballistics).  The  air 
cavity  effect  was  deemed  to  be  minimal  due  to  the  low 
inject  velocities  used. 

C.  Data  Analysis 

The  3-D  data  provided  by  each  camera  was  first  used  to 
produce  raw  two  2-D  plots  of  the  mine's  trajectory.  Next,  2- 
D  data  from  both  cameras  was  then  fused  to  produce  a  3-D 
history.  This  3-D  history  was  then  made  non-dimensional 
in  order  to  generalize  the  results.  The  non-dimensional  data 
was  used  to  generate  impact  scatter  plots  and  was  also  used 
in  multiple  linear  regression  calculations. 

6.  EXPERIMENTAL  RESULTS 

A .  T rajectory  Patterns 

After  analyzing  the  3-D  data  set,  seven  trajectory  patterns 
were  found.  The  plots  on  the  (y,  z)  plane  were  chosen  for 
trajectory  analysis,  as  this  plane  was  parallel  to  the 
direction  of  the  mine  drop.  The  generalized  trajectory 
patterns  are  described  in  Table  3,  and  Figures  7-10.  The 
water  phase  trajectory  a  mine  experiences  ultimately 
determines  the  impact  orientation.  In  MIDEX,  the 
categorizing  of  trajectories  into  general  patterns  served  two 
purposes.  Observed  trajectories  were  found  to  be  most 
sensitive  to  %  -value,  drop  angle  and  LID  ratio.  As  COM 
distance  ( X  -value)  increased  from  GC  the  mine  tended  to 
follow  a  straight  pattern.  As  COM  was  moved  closer  to  the 
GC  (decreasing  x  -value)  the  mine's  trajectory  tended 
towards  being  more  parallel  with  the  pool's  bottom.  At 
steep  drop  angles,  the  mine  experienced  little  lateral 
movement  and  tended  towards  a  straight  pattern. 
Additionally,  as  LID  ratio  decreased  more  complex 
trajectory  patterns  developed.  This  included  significant 


oscillation  about  the  vertical  axis  and  increased  lateral 
movement. 

Table  3.  Description  of  relative  coordinate  based  trajectory 
patterns. 


Trajectory  Pattern 

Description 

Straight  (or  Slant) 

Little  angular  change, 

Nearly  parallel  to  z-axis, 

y/2  near  (90°  ±  15°) 

Spiral 

Oscillating  (no  rotation) 

Flip 

Only  once  as  X  *  <0 

Flat 

IX,  near  0°  (no  oscillation) 

Seesaw 

y/2  nearO0  (oscillation) 

Combination 

Several  of  the  above  patterns 

Straight  Spiral  Flat 

L»  15cm;  Vl»  1  835rVs  L*  9cm:  VS*  1  553mfc  L=  15cm:  2  18nVs 


Figure  7.  Trajectory  patterns. 


Flip 


Seesaw 


Flip -Flat 


L=  1 5cm:  Vi=  2  874m/s  L_  ,  5^,.  Vl=  ,  379m/s 


L=  9cm:  VS*  2  294m/s 


Y<m>  Y(m)  Y(m) 


Figure  8.  Trajectory  patterns. 
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Figure  9.  Trajectory  patterns  for  x  =  X 2  • 

I  I  COM  Position:  -2  I _ 
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Drop  Anqle:  45° 

Flip-Spiral-Slant 

Flip-Slant 

Flip-Spiral 

Flip-Slant 

Flip-Slant 

Flip-Spiral 

Drop  Anqle:  60° 

Flip-Straiqht 

Flip-Straiqht 

Flip-Spiral-Seesaw 

Straiqht-Flip 

Slant-Flip-Slant 

Flip-Spiral 

Figure  10.  Trajectory  patterns  for  negative  COM  position 
( X-\  ?  X-2  )•  Here,  COM  is  above  GC  when  the  cylinder  enters 
the  surface). 

B.  Impact  Attitude 

The  angle,  y/2+7T  I  2 ,  is  the  impact  attitude  at  the  bottom 

of  the  water  column.  The  mine  burial  is  largely  determined 
from  the  impact  attitude  of  the  mine.  Mines  whose  impact 
attitudes  are  perpendicular  (i/f7=  90°)  to  the  sediment 

interface  will  experience  the  largest  degree  of  impact  burial 
(Taber  1999).  It  is  therefore  important  to  analyze  the 
relationship  between  impact  attitude  and  the  controlled 
parameters,  drop  angle,  Vm,  LID ,  and  % .  The  experiment 
shows  that  both  LID  and  V-m  had  little  influence  on  impact 
attitude.  The  drop  angle  and  x  >  however,  were  the  largest 
determinants  of  impact  attitude. 

From  Figure  11  it  is  apparent  that  there  are  several  peaks 
centered  near  90°,  140°,  and  180°.  Further  analysis  reveals 
that  these  peaks  correspond  to  the  COM  positions  0,  1  and 

2  (corresponding  to  Xq  ,  X\  •>  X 2  X  respectively.  COM 

positions  -1  and  -2  (corresponding  to  X-\ ,  X-2)  followed 

the  same  trend  as  their  positive  counterparts. 

Although  drop  angle  was  not  the  most  influential 
parameter,  variations  did  induce  changes  in  impact 
orientation.  As  drop  angle  increased,  the  likelihood  of  any 
lateral  movement  decreased.  This  allowed  for  impact  angles 
that  were  more  vertically  orientated.  This  is  primarily  due 


to  the  fact  that  the  vertical  components  of  velocity  were 
greater  than  those  at  shallow  angles.  Thus,  the  time  to 
bottom  and  time  for  trajectory  alteration  was  less. 

C.  Mine  Tumbling  Not  Observed 

The  current  Navy's  mine  impact  burial  model 
(IMPACT25/28)  was  developed  by  the  Coastal  System 
Station  (Amone  and  Bowen  1980;  Satkowiak  1987); 
subsequent  upgraded  by  the  New  Zealand  Defense 
Establishment  (Hurst  1992)  and  the  Naval  Research 
Laboratory  (NRL).  Some  of  the  major  input  parameters  to 
the  model  are  environment  (sedimentation,  shear  strength, 
water  depth),  mine  characteristics  (shape,  center  of  mass, 
weight,  and  mine  deployment  parameters),  deployment 
platform  (ship,  aircraft,  submarine),  speed  of  platform, 
angle  of  mine  upon  entering  water,  rotational  velocity  at 
time  of  deployment  and  others.  The  theoretical  base  of 
IMPACT25  is  the  momentum  equations  (2a)-(2c).  The 
model  assumes  that  COM  coincides  with  GC  (x  =  0).  Chu 
et  al.  (2001)  reported  the  discrepancy  between  observed 
and  model  predicted  (by  IMPACT25/28)  mine  burial  depth 
from  a  mine  impact  burial  experiment  performed  near  the 
Monterey  beach,  California  in  May  2000.  The  model 
describes  two  trajectory  patterns  (Figure  12):  (a)  without 
any  orientation  change,  (b)  with  a  constant  tumbling.  For 
the  second  pattern,  user  should  input  the  tumbling  rate  (Chu 
et  al.  1998;  2000a, b).  In  MIDEX,  mine  tumbling  was  never 
observed  even  for  x  =  0  (Figure  13). 


Figure  11.  Relationship  between  COM  position  and  impact 
attitude. 
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Figure  12.  Two  trajectory  patterns  predicted  using  the  Navy’s 
Mine  Impact  Burial  Model  (IMPACT25/28):  (a)  no  tumbling, 
and  (b)  constant  tumbling.  The  model  was  established  based 
only  on  the  momentum  equations  (2a)-(2c). 
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Figure  13.  Trajectory  patterns  for  the  COM-O  position  (i.e., 
COM  coincidence  with  GC). 


7.  STATISTICAL  PREDICTION  MODEL 

Multivariate  linear  regression  model  was  established 
from  the  MIDEX  data  (712  drops)  to  establish  relationships 

between  the  input  non-dimensional  parameters;  yf™  ,  LID , 
Fin,  and  x >  and  the  output  variables  (temporally  varying) 
such  as  position  (xm,  ym,  zm),  velocity  ( u ,  v,  w)  and  attitude 
y/2  at  time  t.  Let  Y  represent  the  output  variables.  The 
regression  equation  is  given  by 

r  (0  =  A  (0  +  A  m:  +P2{t)L/D  +  P3  (OF  +  A  (0Z-  (9) 

Figure  14  shows  the  temporally  varying  regression 
coefficients  for  u ,  v,  w,  and  y/2.  For  the  mine's  attitude 

( y/2 )  the  coefficient  /?(0  *s  much  larger  than  the  other 

coefficients,  which  indicates  that  the  mine's  orientation 
(versus  the  vertical  direction)  is  mainly  determined  by  the 
location  of  COM. 

To  show  the  dependence  of  impact  of  mine  on  the 
bottom  (i.e.,  horizontal  location  relative  to  the  mine's 
surface  entry  point,  orientation,  velocity),  the  multivariate 
regression  between  the  input  non-dimensional  parameters; 

y/ 2”  ,  LID ,  Fin,  and  x  5  and  the  final  state  (i.e.,  impact  on  the 
bottom)  variables  such  as  the  horizontal  position  of  COM 
(xm,  ym),  the  velocity  of  COM  ( u ,  v,  w)  and  the  attitude  y/^ . 

Let  Z  represent  the  output  variables.  The  regression 
equation  is  given  by 

Z  =  ao  +  ay"  +  aL  I D  +  aVn  +  x,  (10) 

with  the  regression  coefficients  (ao,a,a2,a,a4)  for  the 
output  parameters  given  by  Table  4. 


Table  4.  Regression  coefficients  of  Eq.(10). 


Xm 

Tm 

v2 

u 

V 

W 

a0 

-0.07 

-0.05 

103.0 

0.0040 

-0.014 

-0.948 

0.119 

-0.83 

-13.4 

-0.008 

-0.0106 

-0.108 

a2 

-0.47 

-0.08 

-0.501 

-0.110 

0.0005 

0.0295 

0.037 

0.062 

1.045 

0.0025 

0.0011 

-0.221 

a4 

0.237 

0.433 

472.0 

-0.009 

0.0537 

-1.25 

8.  CONCLUSIONS 

(1)  Moment  of  a  falling  mine  in  water  column  is  a  highly 
nonlinear  process,  which  should  be  described  by  both  the 
momentum  and  moment  of  momentum  equations.  If  the 
moment  of  momentum  equations  are  absent,  the  motion  of 
falling  mine  cannot  be  completely  simulated.  The  new  mine 
impact  burial  prediction  model  should  include  the  moment 
of  momentum  equations. 

(2)  Six  different  trajectory  patterns  (straight,  spiral,  flip, 
flat,  seesaw,  combination)  were  detected  from  MIDEX.  No 
tumbling  of  mine  was  observed.  The  flip  of  mine  occurs 
only  once  for  negative  x  -values  (i.e.,  COM  is  above  GC  as 
the  mine  enters  the  water  surface).  The  flat  pattern  occurs 
usually  for  x  =  9  (i.e.,  COM  coincides  with  GC).  The 
transition  between  patterns  depends  on  the  initial  conditions 
(drop  angle  yf"  and  initial  velocity  Fin)  and  the  internal 
structure  of  mine  (such  as  LID  ratio,  ^-value,  etc.).  The 
dynamics  of  trajectory  pattern  formation  and  transition  is 
very  complicated.  It  involves  stability,  nonlinear  dynamics, 
and  fluid-body  interaction. 

(3)  MIDEX  shows  that  both  LID  and  Fin  had  little 
influence  on  mine's  impact  attitude  on  the  bottom.  The  drop 

angle  yf™  and  x ,  however,  were  the  determinants  of 
impact  attitude.  For  x  =0,  the  mine  was  almost  parallel  to 
the  bottom.  For  X-2  and  Z2  cases,  the  mine  is  almost 
vertical  to  the  bottom. 

(4)  MIDEX  provided  nondimensional  data  of  position 
(xm,  ym,  zm),  velocity  (1/,  v,  w)  and  attitude  y/2  for  each 

input  including  data  yf"  ,  LID  ratio,  Fin,  and  x  •  The  data 
can  be  used  for  model  development  and  validation. 

(5)  The  observed  trajectories  were  far  more  complex 
than  those  theorized  by  using  only  the  momentum 
equations.  Simply  assigning  a  rotation  rate  into  the  model 
will  not  simulate  the  movement  of  falling  mine  in  the  water 
column.  At  a  minimum,  updates  to  the  IMPACT  25  model 
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should  include  the  more  realistic  moment  of  momentum 
equations. 

(6)  Further  research  on  mine  hydrodynamics  is  needed. 
The  research  needs  to  expand  beyond  the  simple  cylindrical 
shaped  mine  to  those  that  are  irregularly  shaped  (Rockan 
and  Manta  types).  Additionally,  the  utilization  of  scaled 
down  versions  should  be  explored.  A  smaller  mine  that  can 
be  modeled  as  accurately  as  its  real  counterpart  will  save 
time,  money  and  will  require  less  logistical  support. 
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